Fermi Arcs in the Superconducting Clustered State for Underdoped Cuprates 
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The one-particle spectral function of a state formed by superconducting (SC) clusters is studied 
via Monte Carlo techniques. The clusters have similar SC amplitudes but randomly distributed 
phases. This state is stabilized by the competition with antiferromagnetism, after quenched disorder 
is introduced. Fermi arcs between the critical temperature Tc and the cluster formation temperature 
scale T* are observed, similarly as in the pseudogap state of the cuprates. The arcs originate at 
metallic regions in between the neighboring clusters that present large SC phase differences. 



Introduction. Research in hole-doped high tempera- 
ture superconductors is currently mainly focusing on the 
pseudogap (PG) phase that exists above in the under- 
doped regime. This exotic phase acts as the normal state 
to superconductivity in a broad range of hole densities, 
and it may contain the solution to the puzzling properties 
of these compounds. Two groups of experiments have re- 
cently provided important microscopic information about 
the PG phase: (i) Using anglc-rcsolved photoemission 
(ARPES) techniques [ij, the PG Fermi surface was found 
to be composed of disconnected segments, widely known 
as "Fermi arcs" that are centered at the nodal (N) 
points k=(±7r/2, ±7r/2), in the usual two-dimensional 
(2D) square-lattice notation. These arcs are caused by 
the low-energy PG, which presents a k-dependence re- 
sembling a d-wave SC gap close to the anti-nodal (AN) 
points (0, ±7r),(±7r, 0) but it has gapless excitations in 
a finite momentum range near the nodes [1, Q . The AN 
gaps remain approximately constant in the PG phase, in 
contrast with mean- field theories ft'^) Using scan- 

ning tunneling microscopy (STM) techniques, remarkable 
results were recently reported [5|. At temperatures well 
above Tc, the B2212 local density-of-state (LDOS) still 
closely resembles that of a c?-wave superconductor. E.g., 
at optimal doping and r=120K, a d-wave gap is observed 
in clusters that appear to be randomly distributed 0]. 

Theoretically, the PG state is believed to be either 
(i) a. precursor of the SC state, with phase fluctuations 
destroying superconductivity in an homogeneous state 
made out of small Cooper pairs formed by a strong at- 
traction Q, or (ii) caused by other competing orders. 
However, a third possibility was recently proposed Q. 
Via Landau- Ginzburg (LG) calculations, a state with 
nano-scale SC and antiferromagnetic (AF) clusters is sta- 
bilized when quenched disorder (caused, e.g., by chemi- 
cal doping) influences on an otherwise homogeneous state 
with local coexistence of the SC and AF order parame- 
ters. The SC amplitude is robust and approximately the 
same in all the SC clusters. However, the phase of the 
SC order parameter, while uniform within each SC clus- 
ter, randomly changes from cluster to cluster causing the 
overall state to become non-SC. This scenario is concep- 



tually different from (i) and (ii), but contains elements of 
both: the AF competing order is needed for the SC clus- 
ters to form and, once formed, phase fluctuations between 
clusters, occurring even in a weak-coupling BCS regime, 
destroy the SC condensate. This "Superconducting Clus- 
tered State" (sees) can account, at least qualitatively, 
for the STM experiments ja], since both are based on 
a similar inhomogeneous distribution of SC gaps above 
Tc- In this manuscript, it is shown that the SCCS also 
produces a Fermi surface (FS) with Fermi arcs. 

LG simulations. Our calculations start with a LG the- 
ory for the AF vs. c?-wave SC competition, studied with 
Monte Carlo (MC) techniques. Details were extensively 
discussed before 8] and only a brief summary will be 
repeated here. The LG classical Hamiltonian involves a 
complex number Ai=|Ai|e"^' and a real vector Si, rep- 
resenting the SC and AF order parameters at site i of 
a 2D square lattice. The interactions are the standard: 
1. terms with up to 4th powers of the order parameters, 
locally favoring SC and AF; 2. nearest-neighbor (NN) 
couplings that spread the range of the order, emerging 
from gradients in the continuum formulation; and 3. an 
interaction between the order parameters, with strength 
usc|AF- Quenched disorder is also included to represent 
chemically doped cuprates. More specifically, 



H = 



i i 

^PAF(i, a)Si • Si+a +wsc|afX! I^il^l^'l^ 

i,a i 

^psc(i,a)|Ai||Ai+Q,|cos(0i - 

i,a. 



Previous investigations [3] showed that fixing, e.g., r^c — 
-1, rAF=-0.85, usc=WAF=l, -usc|AF=0.7 and varying 
psc and paf along the line paf—^+Psc, produces a 
clean-limit phase diagram with a SC-I-AF region of local 
coexistence Fig. la qualitatively illustrates this clean 
limit case {D=0) [lO|. However, quenched disorder D^O, 



2 



o-o Short Range 
e-oLong Range 




FIG. 1: (Color online) (a) Schematic T vs. doping [x) phase 
diagram 8] for the competition AF vs. SC: _D=0 {D^Q) is 
the clean (dirty) limit. The vertical (green) line is the T- 
range emphasized here: at Tc long-range order develops, and 
at T* SC clusters (short-range order) are formed. In prac- 
tice, the AF to SC transition is reached in Eq.(l)by varying 
psc and paf, keeping paf=1+Psc for simplicity [8[. The dis- 
order in the couplings is correlated with power 0=0.8 3. 
(b) SC correlation (SS) vs. T at short (vectorial distance 
(3,0)) and long (vectorial distance (16,16)) distance, in the 
clean and dirty limits. Results are from MC simulations 
on a 32x32 cluster with periodic boundary conditions, us- 
ing 2,000 (3,000) thermalization (measurement) steps, start- 
ing with a random initial configuration, psc was taken from 
a bimodal distribution with values -1.1 and -0.1 in the dirty 
limit, and uniformly -1.1 in the clean limit. SS is defined as 
■^"^('^^JvITIir l^jll^j+i|cos(0j-(?!)j+i). Some results gath- 
ered on 64x64 lattices show that size effects are not impor- 
tant, (c) Classical SC order parameter A for a typical MC- 
equilibrated configuration at several T's (0.2, l.O^Tc, 1.5 and 
2.0~r*). Color intensity represents |A| and the actual colors 
represent the SC angle (j) at each site (see color wheel). The 
maximum value of |A| for the temperatures studied was ~2.0. 
The bimodal couplings configuration is also shown (d): blue 
indicates regions where SC couplings dominate (psc=— 1-1); 
red where AF couplings dominate (psc =— 0.1). 



for example introduced in Eq.(l) via a random distribu- 
tion of Psc and paf values, reduces both critical temper- 
ature, opening a gap between the competing phases. In 
this glassy region, nano-scale SC and AF clusters coex- 
ist, as sketched in the inset 0]. Two temperature scales 
emerge naturally: upon cooling, first the SC amplitude 



develops at a crossover scale T*, but the SC phases be- 
tween clusters remain random. Reducing T further, co- 
herence among the cluster phases is reached at Tc. The 
presence of two temperature scales is indeed observed 
numerically, as shown in Fig. lb: at D=0, the long and 
short correlations are very similar, while with disorder 
ZJ^O, a substantial difference between T* and Tc (defined 
by the T where short and long range correlations vanish, 
respectively) is clearly observed. MC-equilibrated config- 
urations of Eq.(l) better clarify this issue (Fig. Ic). The 
complex-number SC orders parameters are represented 
by a color and an intensity (see wheel in Fig. Ic). At 
low T, the uniform and intense color indicates a robust 
SC state. As T increases, phase fluctuations appear near 
Tc. In the interesting range Tc<T<T*, both SC clusters 
with random phases and non-SC (white) regions coexist: 
this is the SCCS emphasized here QJ^. Finally, near T* 
or above, few vestiges of SC remain. 

Fermions in the SC+AF background. After equili- 
brated configurations of the SCCS state are gathered 
from the LG/MC classical analysis, fermionic properties 
are obtained by locally coupling itinerant electrons (sim- 
ulating carriers) to the classical order parameters, as pre- 
viously discussed The Hamiltonian is: 



<ij>,'T i 

+ \Y. ^'I^-I' - E ^i(AiaCiTCi+ax + H.c), (2) 



=("iT 



-nii)/2 



where c\a are fermionic operators, s?= 
n-ia is the number operator, and AiQ,=/3| Aile"^' are com- 
plex numbers for the SC order parameter defined now at 
the links (i,i-|-a)(a=unit vector along the x or y direc- 
tions; (3=1 (-1) for a along x (y)). For each {Aid 
{S?} configuration, the fermionic sector is exactly diag- 
onalized via library subroutines and any property, static 
or dynamic, can be easily obtained. At Ji=0, d-wave 
SC is favored since the pairing term involves NN sites, as 
in any standard mean-field Bogoliubov-de-Gennes (BdG) 
approximation . The parameters of relevance are Ji 
and Vi {t is the energy unit), and they carry a site de- 
pendence to easily include quenched disorder 14 1. 

Fermi arcs. Fig. 2 contains our most important results, 
(a) shows the one-particle spectral function ^(k, cu) along 
a straight line from the N to the AN points. At low 
T<Tc, the d-wave SC gap is clearly visible (higher |a;| 
peaks are related with a nodeless AF gap, and they do not 
affect lower energy features). However, as T is raised first 
across Tc, then through the intermediate SCCS phase- 
fluctuating regime proposed for the PC state, and finally 
to above T* , clearly the gaps disappear forming segments 
(arcs) starting at the node, with a length that grows with 
increasing T. In (b), the FS's are shown: four nodes at 
low T become arcs at higher T, they eventually merge, 
and form a closed FS at the highest T. (c) contains an 
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FIG. 2: (Color online) (a) A(k,uj) for equally-spaced k's along 
the direction from {n/2,n/2) (bottom) to (0, vr) (top). The 
classical LG configurations used are obtained as described in 
Fig. 1 but for the case psc=— 0.8 (0) in the SC (AF) regions, 
using paf = 1+Psc, on a 32x32 lattice. The T's are, from 
left to right, T^O.l, 0.4, 1.0 and 2.0. In this case, rc~0.4, 
and T*~ 1.2-1.5. (b) A(k,ui= 0) in the kx-ky plane for the 
same parameters and temperatures as in (a). Results shown 
are those within a window Ae=0.1 from the Fermi level, (c) 
SC gap (distance between peaks) vs. momentum for the same 
parameters as in (a) but with psc^ — 1.1 or —0.1 (bimodal dis- 
tribution), at the T's indicated, (d) A(k,aj= 0) from (0,7r/2) 
to (7r/2,7r), at T=1.0 and parameters as in (a). 



example of gaps vs. k along the N-AN line, showing the 
arc formation, and the stability of the AN gap even when 
T is varied over a wide range. Here, as in [3l, 141], a given 
k is said to have a gap if two peaks are found in A(k, lo). 
In Fig. 3a, the length of the Fermi arc is shown vs. T, 
and an approximate linear relation is found. All these 
results are in good agreement with ARPES [3, 4]. 

An important observation is that the arcs are not 
merely caused by the broadening of the peaks by dis- 
order. To understand this, consider now the direction 
perpendicular to the N-AN line. In Fig. 2d, A{k,uj) 
from (0, 7r/2) to (7r/2, tt) is shown. This crosses the N-AN 
line at (7r/4, 37r/4). The figure shows that in the range 
Tc<T<T* a metallic dispersion, close to non-interacting 
electrons, is observed. This is totally different from the 
low-T results that show a BdG quasiparticle dispersion 
and a gap (not shown). Concomitant with this behav- 
ior, by monitoring the LDOS in our simulations, we have 
noticed the existence of gapless metallic patches coexist- 
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FIG. 3; (Color online) (a) Length of the Fermi arc (as a % of 
the maximum length) vs. T for the case described in Fig. [2^. 
The Fermi arc was defined to exist at a certain momentum k 
when its intensity was within 35% of the maximum intensity. 
Other definitions lead to a similarly linear relation but with 
different T—fO limits. Note that the energy cutoff Ae=0.1 
and the finite momentum resolution, due to the lattice's fi- 
nite size, prevent the size of the nodes from being exactly 
zero in the low-T state. If a shift downwards by 20% (the 
value at T=0) is carried out, an even better agreement with 
[3, 01 is obtained, (b) Schematic representation of the toy 
model configuration (see text), (pi and (j)2 refer to the SC 
order parameter phases in 4x4 squares, (c) LDOS for the 
example shown in (b) with </>i=0 and (j>2='T- Red (blue) solid 
(dashed) lines correspond to a site at the center (border) of 
the 4x4 square. The parameters used are the same at each 
site: |A|=1, V=0.25, and J=2. (d) A(k,uj) for k along the 
direction (7r/2,7r/2) to (0, vr) for the case shown in (b) with 
4>i=4>2 ~ (i.e. perfect d-wave superconductor), (e) Same as 
(d), but with 01 and (j)2 randomly chosen between and vr. 
(f) j4(k, (j;= 0) in the kx-ky plane for case (e) (Ae=0.1). 



ing with the AF and SC clusters. These metallic regions 
appear in "fragile" zones of the disordered configuration, 
such as in long and thin areas of one phase penetrating 
into the other, where none of the two orders prevails. 
Thus, the SCCS actually involves three ingredients: SC, 
AF, and metallic areas. Our numerical simulations sug- 
gest that the metallic areas and Fermi arcs are related. 

Toy models. To better understand the Fermi arc 
formation in the SCCS, simplified models were ana- 
lyzed. Consider a 2D square lattice regularly divided 
into smaller 4x4 squares (Fig. all with the same SC 
amplitude but different phases, and without AF. To sim- 
ulate T effects, frozen configurations {(/>} of the SC phases 
were studied. In Fig. 3d, the uniform case </'i=<^2=0 is 
shown: as expected, a clear d-wave gap exists along the 
N-AN hue. This mimics the regime T<Tc in Fig. [Hi. To 
simulate Tc<T<T* , consider now a random distribution 
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of phases. The result (Fig. ^) still shows a clear gap in 
the AN point, but near the nodal point now a finite set of 
momenta do not present a gap anymore, thus generating 
a Fermi arc (see also the FS in Fig.[3f)). Studying A(k, uj) 
for a variety of 02 's, at fixed <pi—0, we observed that an- 
gles such that cos(02)<O are those that most contribute 
to the arcs. For, e.g., 4'2=t^ metallic portions were iden- 
tified at the lines separating two domains (Fig. 3c) [Tsj l. 

















1^ 


^2 



(a) 






-4 


-2 2 


4 




CO 


















-4 -2 2 4 







Conclusions. A state based on SC clusters with ran- 
dom phases was here used to calculate the ARPES re- 
sponse of the PG state. Fermi arcs were found in a T 
range between Tc and the T* where the clusters start 
forming upon cooling. Toy models illustrate and sim- 
plify these results. For further progress, our effort must 
be supplemented by better analytic control of the SCCS 
and by numerical studies using microscopic models. The 
extension of these results to s-wave superconducting films 
should also be pursued [l^ . In addition, the STM exper- 
iments @ not only showed a clustered stated in a broad 
temperature range, but also unveiled a regime immedi- 
ately above where 100% of the sample had the d-wave 
gap. Adapting the SCCS to also accommodate this fea- 
ture will require increasing the attraction V [13] and this 
direction will be the focus of future investigations. 
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FIG. 4: (Color online) (a) Schematic representation of a 
Josephson-junction-like structure. Top: two SC regions with 
order parameters Ai and A2, separated by a A=0 area of 
width w (grey). Middle: expected gap interpolation due 
to proximity effect, for Ai=A2. Bottom: same as middle, 
but for Ai=-A2. (b) LDOS on a 32x32 lattice containing 
4 equally-spaced 12x12 SC clusters separated by ^1=4. The 
order parameters are staggered between Ai=l and A2=-l, 
with y=0.25 and without AF (J=0). Shown are results from 
the center of one SC region to a nearest-neighbor, (c) Den- 
sity of states of a perfect superconductor (blue solid) and of 
a perfect metal (black dashed) on a 32x32 lattice. Note the 
Van Hove singularity in the 2D metal, (d) A{\i,uj = 0) in 
the kx-ky plane for Ai = — A2=l and w=2, with a setup as 
described in (b) but using 14x14 SC clusters. 

Relevance of large phase differences. The metallic re- 
gions seem caused, at least in part, by large phase differ- 
ences between neighboring SC clusters. Consider Fig.|4^: 
here two SC clusters are shown, with order parameters 
AiT^O and A27^0, separated by an intermediate thin re- 
gion where A=0. If Ai=A2, this intermediate region will 
develop a gap by proximity effect (Fig. 2^), and the FS 
has still 4 nodes, as confirmed by an explicit BdG cal- 
culation (not shown). However, if Ai=-A2 a qualitative 
difference occurs: now the interpolation in the interme- 
diate region necessarily requires the existence of a zero, 
where the SC order parameter must vanish even for a 
thin A=0 layer (Fig. gla,). A BdG study confirms that 
the LDOS in the middle between two "anti-parallel" SC 
clusters is almost identical to the LDOS of a metallic 
state (see Figs. 4b, c). The FS for Ai=-A2 has Fermi 
arcs (Fig. 4d). Thus, large differences induce metallic 
regions, and those appear to cause the Fermi arcs. 
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